fix(.First)
options(editor="notepad")
fix(.First)
fix(.First)
.First()
search()
fix(.First)
?PCAbipl
search()
.First()
search()
?CVAbipl
fix(matmult)
package.skeleton(name = "MatMFor", list= "matm", path =  
 "G:\\Workshop 25-04-2012", namespace = TRUE)
package.skeleton(name = "MatMFor", list= "matmult", path =  
 "G:\\Workshop 25-04-2012", namespace = TRUE)
package.skeleton(name = "MatMFor", list= "matmult", path =  
 "G:\\Workshop 25-04-2012\\MatMFor", namespace = TRUE)
objects()
dump("matmult", file = "G:\\Workshop 25-04-2012\\matmult.rfun")
fix(.First)
.First()
matmult
objects()
fix(matmult.Fortran)
A <- matrix(rnorm(100) nrow=20,ncol=5)
A <- matrix(rnorm(100), nrow=20,ncol=5)
B <- matrix(matmult.Fortran<-function (A,B) 
{if(ncol(A) != nrow(B)) stop("A and B not conformable with respect to matrix multiplication \n")
n <- nrow(A)
q <- ncol(B)
p <- ncol(A)
Cmat <- matrix(0, nrow=n, ncol=q)
storage.mode(A) <- "double"
storage.mode(B) <- "double"
storage.mode(Cmat) <- "double"
value <- .Fortran("matm", as.integer(n), as.integer(p),as.integer(q),A,B,matprod=Cmat)
value$matprod
}
B <- matrix(rnorm(25) nrow=5,ncol=5)
B <- matrix(rnorm(25), nrow=5,ncol=5)
matmult.Fortran(A,B)
matm(A,B)
matmult(A,B)
sys.time()
system.time()
?system.time
system.time(matmult(A,B))
B <- matrix(rnorm(2500), nrow=50,ncol=50)
A <- matrix(rnorm(10000), nrow=200,ncol=50)
system.time(matmult(A,B))
system.time(matmult.Fort(A,B))
system.time(matmult.Fortran(A,B))
A <- matrix(rnorm(100000), nrow=2000,ncol=50)
system.time(matmult.Fort(A,B))
system.time(matmult.Fortran(A,B))
system.time(matmult.Fortran(A,B))
system.time(matmult.Fortran(A,B))
system.time(matmult.Fortran(A,B))
system.time(matmult.Fortran(A,B))
system.time(matmult.Fortran(A,B))
system.time(matmult.Fortran(A,B))
system.time(matmult(A,B))
system.time(matmult(A,B))
system.time(matmult(A,B))
?system.time
proc.time()
proc.time()
Fortran.example <- function(n,p)
{time.begin <- proc.time()[3]
A <- matrix(rnorm(100000), nrow=n,ncol=p)
A <- matrix(rnorm(100000), nrow=n,ncol=p)
}
fix(Fortran.example)
Fortran.example()
fix(Fortran.example)
Fortran.example()
Fortran.example(10,5)
Fortran.example()
fix(Fortran.example)
Fortran.example(10,5)
fix(Fortran.example)
Fortran.example(10,5)
Fortran.example(10000,5000)
Fortran.example(100,50)
Fortran.example(1000,500)
utils:::menuInstallLocal()
utils:::menuInstallLocal()
.First
.First
fix(.First)
options(editor="notepad")
fix(.First)
Fig.2.12()
Fig.2.5()
Fig.2.13()
Fig.2.14()
Fig.2.15()
Fig.2.18()
Fig.2.19()
Fig.2.28()
Fig.2.35()
Fig.2.38()
Fig.2.40()
objects()
Fortran.example()
fix(Fortran.example)
Fortran.example()
Fig.2.22()
args(vectorsum.interp)
vectorsum.interp(4)
Fig.4.18()
Fig.4.18
pCAbipl(CVAbipl.density(CopperFroth.data[, 1:8])
PCAbipl(CVAbipl.density(CopperFroth.data[, 1:8])
PCAbipl.density(CopperFroth.data[, 1:8])
PCAbipl.density(CopperFroth.data[, 1:8],scale=T)
CVA.predictivities()
args(CVApredictivities)
args(CVA.predictivities)
Copperfroth.data
objects(4,pat="data")
CopperFroth.data[1:5,]
CVA.predictivities(CopperFroth.data[,-9],indmat(CopperFroth.data[,9])
)
graphics.off()
Fig.4.18()
Fig.4.18
Fig.4.18
Fig.4.18
Fig.4.19()
Fig.4.20()
4.22()
Fig.4.22()
Fig.4.24()
Fig.4.24()
Fig.4.24()
Fig.4.25()
Fig.4.26()
Fig.4.26
